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ABSTRACT 

We present results from Chandra and XMM-Newton X-ray observations of 
NGC4214, a nearby dwarf starburst galaxy containing several young regions of very 
active star-formation. Starburst regions are known to be associated with diffuse X-ray 
emission, and in this case the X-ray emission from the galaxy shows an interesting 
morphological structure within the galaxy, clearly associated with the central regions 
of active star-formation. Of the two main regions of star formation in this galaxy, 
X-ray emission associated with the older is identified whereas little is detected from 
the younger, providing an insight into the evolutionary process of the formation of 
superbubbles around young stellar clusters. The spectra of the diffuse emission from 
the galaxy can be fitted with a two temperature component thermal model with 
kT = 0.14 keV and 0.52 keV, and analysis of this emission suggests that NGC4214 
will suffer a blow-out in the future. 

The point source population of the galaxy has an X-ray luminosity function with 
a slope of —0.76. This result, together with those for other dwarf starburst galaxies 
(NGC4449 and NGC 5253), was added to a sample of luminosity functions for spiral 
and starburst galaxies. The slope of the luminosity function of dwarf starbursts is seen 
to be similar to that of their larger counterparts and clearly flatter than those seen in 
spirals. Further comparisons between the luminosity functions of starbursts and spiral 
galaxies are also made. 

Key words: ISM: jets and outflows - galaxies: individual: NGC 4214 - galaxies: 
starburst - X-rays: galaxies. 



1 INTRODUCTION 

Hierarchical models of structure formation in the Universe 
suggest tha t dwarf galaxies were some of t he first objects to 
be formed (iNavarro. Frenk fc Whitelll995l) . and hence were 
amongst the earliest sites of star formation. If dwarf galaxies 
in the local Universe are analogues of those formed at high 
redshift then studying the structure and evolution of local 
dwarf galaxies will be useful in understanding star formation 
at high redshifts. 

When starburst galaxies undergo intense periods of star 
formation many supernovae occur in a relatively small re- 
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gion of the galaxy. The ejecta from these combine with the 
winds and ejected material from the massive stars in OB 
associations and Super Star Clusters (SSCs) in the sur- 
rounding area to form a superbubble around the cluster 
which expands, sweeping up surrounding medium to pro- 
duce large volumes of hot, shock e d gas with T ~ 10 8 K. 
The simulations of ISuchkov et all <ll994l) show that X-ray 
emission from starburst galaxies is most likely to be due to 
the shocks formed from the interaction between the star- 
burst outflow and the surrounding medium, rather than 
from the material within the winds. The bubbles of hot gas 
expand along the steepest density gradient - usually the mi- 
nor axis of the galaxy, and will eventually expand out of the 
galaxy's gravitational potential well. Indeed, Ha images of 
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starburst galax ies show filaments extending along the m i- 
nor axis jMarlowe et alJll995l iLehnert fc Heckmanlll996bl) . 
The low mass and hence small gravitational potential well of 
dwarf galaxies means that it more likely than in other galax- 
ies for material to escape into the surrounding intergalactic 
medium and form galactic superwinds. In some models these 
starburst-driven winds result in enough interstellar medium 
(ISM) being ejected from the galaxy via galactic super- 
winds to prevent any further star formation from being able 
to occu r foe Young fe Heckmanll994T l. lLehnert fc Heckmanl 
(1996a) find that superwinds are most commonly formed in 
galaxies which have extreme IR properties, such as high IR 
luminosities and warm far-IR flux. Dwarf starburst galaxies 
can therefore play an important role in enriching their sur- 
roundings with heavy elements produced during sta rburst 
episodes foerrara fc Tolstovll200a iDekel fc Silkll 19861) . 

NGC 4214 is a nearby dwarf starburst galaxy with an 
almost face-on orientation. Using the tip of the red gi- 
ant branch method the distance to NGC 4214 has b een 
estimated bv lMafz-Apellaniz. Cieza fc MacKentvl i2002T) as 
D = 2.94Mpc, and this is the value used throughout this 
paper. We note that a d istance of 2.7Mpc was derived by 
iDrozdovskv et alJ (120021) using the same method. A few 
x 10 r years ago the galaxy underwent a burst of star for- 
mation during whic h around 5% of the stellar mass of the 
galaxy was formed (iHuchra et alj ll983l . Colour magnitude 
diagrams for the field stars outside the main star-forming re- 
gions show evidence for both old and young populations, and 
are typical for dwarf galaxies. The diagrams are dominated 
by the "red tangle" - an area of the diagram that contains 
RGB, AGB and blue-loop stars - which implies that there 
is an underlying population of old stars iDrozdovskv et alJ 
2002). The basic properties of NGC 4214 are given in Ta- 
ble □ 

NGC 4214 has been studied at many wavelengths, each 
of which provide an insight into its structure and composi- 
tion. These observations appear to be show ambiguous re- 
sults when considering whether NGC 4214 has undergone a 
recent interaction with a nearby galaxy. / band observations 
of NG C 4214 show a smooth, symmetric disk dFanelli et alJ 
119971) . The outer contours of the I band emission are regu- 
lar, which suggest that there is no evidence that the galaxy 
has undergone a recent merger. NGC 4214 is Hi-rich like 
many other dw arf starbursts. In contrast, observations of 
HI emission by IA11sopp| lll979l) show an extension towards 
a nearby irregular galaxy NGC 4190 suggesting a possible 
tidal interaction between the two galaxies could have trig- 
gered the starburst in NGC 4214. Observations in the far 
ultra-violet (FUV) shows the presence of knots of emission 
at the centre of th e galaxy embedded in diffuse emission 
foanelli et al.lll997l) . These knots denote the presence of 
OB associations which have also been observed at optical 
wavelengths (e.g. IHuchra et aljH983l iLeitherer et al]fl996t 
iMacKentv et alj|2000l) ~ 

More recent observations of NGC 4214 using the Hubble 
Space Tel escope (HST) have show n evidence for a two-stage 
starburst iMacKentv et alJl2000t) . producing SSCs and OB 
associations. These observations show that the galaxy has 
a complicated morphology and is dominated by two large 
HII complexes often referred to as the North West (NW) 
and South East (SE) complexes (these structures are also 
referred to as NGC 4214-1 and NGC 4214-11 in the literature 



Table 1. Properties of NGC 4214 



Parameter 


Value 


Ref. 


a [2000] 


12 h 15 m 39.26 s 




8 [2000] 


+36°19'33.4" 




Distance 


2.94Mpc 


1 


1 arcsec 


855pc 




D 25 


8.91' x 5.62' 


2 


Position angle 


42.2° 


2 


Metallicity 


~ 0.25 


3 


m B 


10.20 mag 


4 




17.9 Jy 


5 


SlOOM"t 


29.0 Jy 


5 


Lfir 


g.lxlO^ergs- 1 


6 


SFR 


0.04 Mq yr _1 


7 



References: (1) IMafz-Apellaniz et al.l <2002l): (2) LEDA 
datab ase (available at http://lcda.univ-lyonl.frl; (3) 
iKobulnickv fc Skillmai] l ll996ft : (4) NED (NASA Extragalactic 
Database, availab le at http://nedwww.ipac.caltech.edul; (5) 
ISoifer et af] <1989ft : (6) see Sectio n IO TableElfor more details; 
(7) using the Lpm relation from iKennicut J Jl998t) . 



as well). In addition many knots of Ha emission are ob- 
served along with structurally varying diffuse gas surround- 
ing the main regio ns of emission. Recent CO observations 
J Walter et all2 001) show three main sites of molecular emis- 
sion, two of which correspond to the NW and SE complexes, 
with an additional site at which star formation has yet to 
commence. 

In optical appearance NGC 4214 is simil ar to other 
dwarf irregular galaxies, however IMafz-Apellaniz et alJ 
(1999) suggest that it has an unusually thin galactic disk 
with a thickness of around 200pc; outflows from the galaxy 
extending this by a further 200pc or so. Flatter galaxies pref- 
erentially undergo blow-out which is when supernovae explo- 
sions blow through the disk of the galaxy, forming a way for 
gas from subsequent explosions to escape without having 
much affect on the surrounding gas foe Young fc Heckmanl 
11994 iMac Low fc Ferraralll9991) . 

The launch of the Chandra Observatory has revolu- 
tionised X-ray astronomy, providing an unprecedented spa- 
tial resolution of 0.5". This resolution enables the identifica- 
tion of point sources down to fainter luminosities than those 
seen with the previous X-ray telescopes such as ROSAT and 
Einstein. It is now possible to study the diffuse X-ray emis- 
sion from galaxies without severe contamination from unre- 
solved point sources, allowing a more accurate definition of 
its properties. The XMM-Newton telescope launched in 1999 
combines both reasonable spatial resolution with very high 
sensitivity. This allows the temperature and the elemental 
abundances of hot, diffuse gas to be studied in great detail. 
The combination of data from both instruments promises a 
much deeper understanding of starburst galaxies. 

Prior to the launch of Chandra and XMM-Newton 
X-ray data from several dwarf starburst galaxies had 
been obtained from ROSAT observations (examples of 
R OSAT data on dwarf starb urst gal axies can be found 
in IStevens fc Str icklan 1 lll998al) and IStevens fc Strickland! 
l|l998bl) '). However. ROSAT lacked both the sensitivity and 
resolution to be able to disentangle the point source emis- 
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sion from X-ray binaries from the diffuse emission, which 
severely limits the interpretation of the data. 

This analysis of NGC4214 adds to a small sample 
of published results from Chandra and XMM-Newton ob- 
servations of dwarf starbursts. One significant example 
is the analysis of a Chandra observ ation of NGC 1569 
jMartin. Kobulnickv fc Heckmanl 12002). With the resolu- 
tion of Chandra it is possible to study not just the global 
properties of the X-ray emiss ion but also t h e pro perties of 
different spatial components. iMartin et al.l J2002T) found a 
strong correlation between the X-ray and Ha emission, and 
found a high a-element to Iron abundance ratio (a/Fe), 
which is consistent with a large amount of Type II supernova 
ejecta contribution iMartin et all2002T) . It is concluded from 
this study that elements can be released from dwarf galaxies 
into the surrounding intergalactic medium. Further compar- 
ison is presented in Section 15.31 In addition to NGC 1569 
published data is available from Chandra observations o f 
the dwarf starbursts NG C 4449 JSummers et alJ l2003aT). 
NGC 3 077 llOtt et alJl2003T) and NGC 5253 iSummers et all 
l2003bl) . 

In Section [5] we describe the Chandra and XMM- 
Newton observations of NGC 4214 and the data reduction 
methods used. In Section |H] the properties of the point 
sources are analysed, and in Section [I] the properties of the 
diffuse X-ray emission are investigated. The morphology of 
this emission is discussed and compared the Ha emission 
from the central regions of the galaxy. In Section |K| the lu- 
minosity function of NGC 4214 is compared to those of spi- 
ral and other starburst galaxies and the properties of the 
diffuse emission and their implications for the formation of 
superbubbles are discussed. Finally in Section |S| the main 
conclusions are summarised. 



2 OBSERVATIONS AND ANALYSIS 
2.1 Chandra Observations 

NGC 4214 was observed using the Advanced CCD Imag- 
ing Spectrometer (ACIS) instrument on board the Chandra 
satellite on 22nd October 2001, using the ACIS-S configu- 
ration. The observation was a single exposure of 26.4ksec, 
although cleaning the data left just over 16.5ksec. The en- 
tire galaxy is included on the back-illuminated ACIS-S3 chip 
hence only this chip was used in the analysis. The data 
reduction was undertaken using Chandra Interactive Anal- 
ysis of Observations (CIAO) version 2.2.1 software analy- 
sis tools . The level 2 event files were filtered to reject all 
events with energies outside the range of 0.3-10.0 keV and 
to include event grades 0,2,3,4 and 6. The background light 
curve was analysed, and all periods of the observation con- 
taining flaring or periods of high background above the 3<r 
level were removed using the lc-clean.sl script provided by 
the Chandra X-ray Center (CXC). 

The point sources were located using the wavdetect pro- 
gramme, and spectral analysis was carried out using XSPEC 
(version 11.2.0). To correct for the low energy quantum effi- 
ciency (QE) degradation caused by molecular contamination 

1 see Chandra data analysis methods, published by the Chandra 
X-ray Center on |http : //asc . harvard . edu/ciao/ 1 



of the ACIS optical blocking filters JPlucinskv et al]l2003h 
the CXC supplied corrarf code was used to correct the ef- 
fective area file. This is dependent on the number of days 
between the launch date of Chandra and the observation, in 
this case day 818. 

Fig. shows the data from the ACIS-S3 chip of the 
NGC 4214 field. The D 25 ellipse of NGC 4214 has semi-major 
and semi-minor axes of 4.5' and 2.8' respectively, corre- 
sponding to an area of ~ 39 square arcminutes. The D25 el- 
lipse is shown in Fig.0 as are the 20 detected point sources. 
The detection and analysis of the point sources are discussed 
in more detail in Section |!|] while the diffuse emission asso- 
ciated with the central regions of NGC 4214 is discussed in 
Section 0] 

Zooming into the central regions of NGC 4214, Fig. 
shows the Chandra X-ray data for the central 2' x 2'. For this 
image the data has been binned by a factor of 2 and then 
adaptively smoothed. This central region contains 3 point 
sources (Srcs 10, 11 and 14), as well as diffuse emission. 
The extent of the NW complex, as seen in Ha emission, is 
depicted by the large circle in Fig. 01 The smaller blue circle 
represents the extent of the SE complex of Ha emission, 
which does not show as much X-ray emission. The central X- 
ray emission, and correlations with Ha emission is discussed 
in detail in Section [4.21 

2.2 XMM-Newton Observations 

The XMM-Newton observation of NGC 4214 was taken on 
November 22nd - 23rd 2001 for a total duration of 14ksec. 
The data was analysed using the XMM-Newton Science 
Analysis System (SAS) version 5.3.3 standard processing 
analysis. The data from the two MOS cameras and the 
PN were processed using the the processing chains emchain 
(MOS) and epchain (PN) to remove both bad pixels and 
bad events, and the data were checked for periods contain- 
ing background flaring. After these periods were removed 
there was around 13ksec of data left for the two MOS cam- 
eras, and 12ksec for the PN. 

When the Chandra and XMM-Newton observations 
were compared it was clear that there was a significant po- 
sitional offset of around 15". To investigate the pointing ac- 
curacy of the observations sources common to both obser- 
vations were correlated to sources found in the APM (Auto- 
matic Plate Measuring) catalogue 2 (including sources found 
on chips other than the ACIS-S3 chip used in this paper). 
Two matches were found, and the positions of the sources 
in the Chandra observation matched those in the APM very 
well, whereas both of the XMM-Newton sources displayed 
and offset of around 15". In addition an archived ROSAT 
observation shows increased X-ray emission corresponding 
with the position of the X-ray emission observed by Chan- 
dra. Throughout this paper the positions obtained from the 
Chandra data are adopted, and the XMM-Newton data has 
been shifted accordingly. 

In Fig. |21 we show a mosaiced image from the 
three XMM-Newton cameras of the central regions of the 
NGC 4214 field. The diffuse emission seen in the Chandra 
data is also visible in this observation, and the emission 

2 Available at http://www.ast.cam.ac.uk/~mike/apmcat/ 
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Figure 1. The NGC4214 field, as seen with Chandra, showing the whole of the ACIS-S3 chip (8.5' X 8.5'), binned up by a factor of 
2. NGC4214 can be seen as a region of diffuse emission at the centre of the image. The large ellipse represents the D25 ellipse of the 
galaxy; the smaller ellipses show the location and extent of the discrete sources as found by wavdetect. 



is clearly not symmetric, with an extension to the North 
East. In the Chandra observation 3 point sources were iden- 
tified in the central regions of NGC 4214, while in the XMM- 
Newton data they are much less clearly resolved. Because of 
the confusion between point sources and diffuse emission in 
the central regions of NGC 4214 we shall concentrate on the 
Chandra data for much of the paper. 



3 POINT SOURCE ANALYSIS 

In this section we describe the detection and investigation 
of the point sources in the NGC 4214 field. The point source 
detection tool wavdetect was used on the Chandra ACIS-S3 
chip, with a detection threshold of 5.41 x 10 -7 (proportional 
to the inverse of the number of pixels in the image) - this 
level means that approximately 1 spurious point source de- 
tection is expected in the field. Using this method, a total 
of 20 point sources were found (as shown in Fig. and 
the positions of these sources are listed in Tabled Of these 
sources a total of 17 fall within the D25 ellipse of the galaxy 
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Figure 2. The X-ray emission from NGC4214 observed with 
Chandra, binned by a factor of 2 and adaptively smoothed. The 
region shown is 2' X 2' . The X-ray emission can be seen to coincide 
with the Ha emission (see text), the NW and SE complexes of 
which are depicted by the two circles. The three central point 
sources, Srcs 10, 11 and 14, are clearly seen. See Section [O] for 
further discussion on the correlation between the X-ray and Ha 
emission. 




Figure 3. The diffuse emission from NGC4214 as observed with 
XMM-Newton. The data is from a mosaiced image of the MOS1, 
MOS2 and PN instruments. The region shown is 2' X 2', the same 
as that of the Chandra image shown in Fig. [5] In the central 
regions of the galaxy at least 2 of the 3 central point sources seen 
with Chandra are apparent, though blended with diffuse emission. 
In addition, the diffuse emission can be seen to extend towards 
the East. 



which is depicted as the large ellipse in Fig.0 Those sources 
outside the D25 ellipse are Srcs 1, 2 and 7. It is worth noting 
that the apparent elongated shape of 1 and 2 are due to the 
shape of the instruments' PSF towards the edge of the chip. 

The point source detection regions generated by wavde- 
tect are output as ellipses which contain 99.7% (3a) of the 
source counts. For each detection a background region was 
defined as an annulus centred on the source with an inner 
radius just larger than the maximum axis of the detection 
ellipse and an outer radius 3 pixels larger than this. Local 
background regions are useful as they allow for any local dif- 
fuse emission on which the point source emission is superim- 
posed to be taken into account. The background subtracted 
counts of the sources are also listed in Table [5] 

For all of the detected point sources, the detected counts 
were extracted for four different energy bands; firstly over 
the whole energy range (0.3-8.0 keV), then in a soft (S) band 
(0.3-0.8 keV) a medium (M) band (0.8-2.0 keV) and finally a 
hard (H) band (2.0-8.0 keV). Background counts were taken 
from the background annulus corresponding to each source 
for each band, and were subtracted from the total. The back- 
ground subtracted values for each band are shown in Tabled 
For two sources, namely Srcs 2 and 8, there is a negative 
number of counts in one of the energy bands due to the 
poor statistics. 

Just 4 sources have > 60 counts; Srcs 2, 10, 11 and 15. 
This is sufficient for a spectral fit to be made, and these 
sources are discussed further in Section ^. 21 For the remain- 
ing 16 sources the Portable, Interactive Multi-Mission Sim- 
ulator, PIMMS 3 , was used to predict a flux assuming the 
exponent of the power law component r=1.8 and a col- 
umn density, Nh, equal to the Galactic value of Nh = 
1.5 x 10 20 cm -2 . 

We find that the luminosity lower-limit for a reliable 
detection in the Chandra observation, assuming that the 
source is at the distance of NGC4214, is Lx ~ 10 36 erg s" 1 . 

Comparing the point sources in the NGC 4214 observa- 
tion with the results from the Chandra Deep Field-South 
observation reported on bv lRosati et alJ i2002f) . we expect 
on average ~ 5 of the sources seen in this observation to 
be background objects and ~ 3 of the sources within the 
NGC 4214 D25 ellipse to be background objects. This sug- 
gests that the majority of sources detected in this observa- 
tion are associated with NGC 4214. 



3.1 Hardness Ratios 

Most of the sources detected in this observation are very 
faint; there are only two sources with > 150 counts and 
just two more with > 60. However, for faint sources it is 
possible to derive some spectral information using hardness 
ratios. Hardness ratios provide information about both the 
absorption and the temperature of point sources. For data 
of sufficient quality it is possible to use a hardness ratio plot 
to distinguish between different classes of sources such as 
super nova remnants, s u perso ft sources, X-ray binaries and 
AGN; IPrestwich et alJ (120031) describe the classification of 
sources using this method. 

The counts in the S, M and H bands were used 



3 Available at |http: / / asc.haivard.edu /toolkit / pimm s.jsp| 
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Table 2. The positions of the 20 point sources detected in the Chandra observation of the NGC4214 field. In addition the number 
of counts (background subtracted) detected for each source in the total energy range and in the soft, medium and hard energy bands 
are shown. Errors are calculated on the counts using Poissonian statistics. For sources with < 60 counts the unabsorbed flux has been 
calculated by assuming an absorbed power law model with an exponent of 1.8 and an absorption equivalent to the Galactic value of 
Njf = 1.5 X 10 20 cm — 2 (see text for full explanation). The fluxes of the 4 sources with > 60 counts have been obtained from individual 
fits to the sources (see Section lij. 21 . 



RA (h m s) Dec (° ' ") net counts Soft band Medium band Hard band Flux 



Source 




J2000 




J2000 




(i error) 
0.3-8.0 keV 


(zb error) 
0.3-0.8 keV 


(i error) 
0.8-2.0 keV 


(zb error) 
2.0-8.0 keV 


erg cm 2 s 1 
0.3-8.0 keV 


ja 


12: 


.ID- 


.11. 


17 


36:20:07. 


.0 


27.0±5.2 
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Q 79. v 10 

y. i o x iu 
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2 a,b 
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1 c 
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l.Oil.O 


2.8±1.7 


6.8±2.7 


3.8 x 10- 


-15 


17 


12: 


15 


:45. 


90 


36:21:37. 


.5 


7.8±2.8 




5.8±2.5 


2. Oil. 4 


2.8 x 10- 


-15 


18 


12: 


15 


:46. 


61 


36:20:07 


.8 


10.0±3.2 


10.0±3.2 






3.6 x 10" 


-15 


19 


12: 


15 


:47. 


46 


36:21:25. 


.4 


11.0±3.3 


4.0±2.0 


4.0±2.0 


3.0±1.7 


4.0 x 10- 


-15 


20 


12: 


15 


: 19. 


24 


36:21:45. 


.7 


11.6±3.5 


l.Oil.O 


1.6±1.5 


9.0±3.0 


4.2 x 10- 


-15 



a sources which fall outside the D25 ellipse of the galaxy. 
b flux obtained from fitted spectrum, see Section l3.2l 



to calculate 2 hardness ratios using ( H— M)/(H+M) and 
(M-S)/(M+S) iSummers et al.1 l2003al) . Fig. g| shows the 
hardness ratio plot. Two sources had no counts in 2 of the 
energy bands, and hence are not included in the plot; neither 
were the two sources for which negative values were found 
in one of the energy bands. 

PIMMS was used to predict colours for point sources 
with both thermal bremsstrahlung and power law mod- 
els, with varying kT and T respectively. In addition, mod- 
els have been calculated with both Galactic absorption, 
(Nh = 1-5 x 10 20 cm -2 ), and with the column density de- 
rived from a spectral fit of the X-ray emission from the entire 
D 25 ellipse of NGC4214 (N H = 1.6x 10 21 cm -2 ). The details 
and further parameters obtained from this fit are discussed 
in Section f5. 21 The theoretical tracks produced by these re- 
sults are also shown in Fig. 2] 

There is no real correlation seen between the two hard- 
ness ratios. There appear to be two distinct groups of 
sources; those with the smallest errors correspond to the 
sources with the largest number of counts and fall in the 
same region of the diagram. These points all fall to the right 
of the track obtained using the higher absorption compo- 
nent implying that Nh = 1.6 x 10 21 cm -2 is a lower limit on 
the value of absorption for each source. None of the sources 
within the D25 ellipse are dominated by counts in the hard 
band, and there are 4 sources that have no counts in the 
hard band at all. This information agrees with the previous 
assertions that the majority of the point sources observed 



are associated with the galaxy and are not objects such as 
background quasars. 



3.2 Individual Point Sources 

Four sources (Srcs 2, 10, 11 and 15) have > 60 counts, and 
these were fitted individually using a power law model with 
an absorption component frozen at the Galactic value of 
Nh = 1.5 x 10 20 cm -2 and with a further absorption compo- 
nent allowed to fit freely. The spectra were grouped so that 
there were a minimum of 5 counts in each bin, with a mini- 
mum of 10 for Src 11. Srcs 2 and 10 do not have sufficient 
counts to allow the power law component to fit freely, and 
hence it was frozen at a value of F = 1.8. In the analysis of 
sources 11 and 15 this parameter was allowed to vary freely. 
The results of the spectral fits are shown in Table |21 

Cross-correlating with the APM catalogue 1 optical 
source was found to correspond to an X-ray detection. It 
corresponds to Src 15, at a separation of 0.22" in RA 
and 0.98" in Dec. The position of the optical source was 
at a[2000] = 12 h 15 m 41.38 s , 5[2000] = 36° 21' 15.7", com- 
pared to a[2000] = 12 h 15 m 41.40 s , £[2000] = 36°21'14.7" for 
Src 15. 

Src 10 corresponds to th e posit ion of a weak Ha source 
observed bv lMacKentv et al.1 feOOOft (denoted as NGC 4214- 
IIIs in th at paper). A source at this position was also ob- 
served by iFanelli et alJ lll997t) in an / band image of the 
galaxy, but not in the FUV. It is likely that this source is 
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Table 3. The components of the models fit to the 4 point sources with > 60 counts. Each was fitted with an absorbed power law model 
with two absorption components; one frozen at the Galactic absorption value of Njj = 1.5 X 10 20 cm — 2 and the other allowed to vary 
to account for any additional local absorption. Column 2 shows the value of the second fitted absorption component and column 3 the 
exponent from the power law component. Column 4 gives the statistics for the fit, and column 5 the absorption corrected luminosity for 
each source. 



Source 


N H (cm" 2 ) 


r 


X 2 /d.o.f. 


Lx (ergs 1 ) 


2 

10 
11 
15 


(8.0±^ (i 2 ) x 10 21 
(1.8±J-|) x 10 21 
(1.8±[] : f ) x 10 21 
(l.0±li) x 10 21 


1.8(frz) 
1.8(frz) 

i.8±g:i 
i-4±8:l 


0.57/10 
1.12/10 
0.8/44 
1.2/24 


(6.0±| |) x 10 37 
(3.6±J :2 ) x 10 37 
(3.0±° : |) x 10 38 
(7.3±f :2 ) x 10 37 



nil = 1.6xl0 a, c 



+ 
•x 



r / 



(M-S)/(M+S) 

Figure 4. The hardness ratios calculated for the point sources 
within the D25 ellipse, omitting 4 sources (see Section l3.1l for ex- 
planation). The energy bands used are: S=0.3-0.8 keV, M=0.8- 
2.0 kcV and H=2.0— 8.0 keV. The connected triangles represent the 
colours predicted by an absorbed thermal bremsstrahlung model 
with kT varying from 0.2 (bottom left) to 2.0 keV in increments 
of 0.2 keV, and from 2.0 to 3.6 keV in increments of 0.4 keV. The 
connected squares represent the colours predicted from an ab- 
sorbed power law model with an exponent T varying from 0.0 
(top right) to 4.0 in increments of 0.4. Models have been calcu- 
lated using the Galactic absorption (Njj = 1.5 X 10 20 cm -2 ), and 
also with the Galactic absorption plus the absorbing column den- 
sity for NGC4214 (N H = 1.6 X 10 21 cm~ 2 , see text for further 
explanation). 



at the dynamical cent re of the galaxy, a view supported by 
the HI observations of iMdntvrel (119981) which place the ro- 
tational centre of the g alaxy close to the position of Src 10 . 

For a comparison, iGrimm. Gilfanov fc Sunvaevl (120031) 
suggest that the number of high mass X-ray binaries present 
in a galaxy is proportional to the star formation rate, and 
use data from a number of galaxies to derive the following 
relation between the number of bright X-ray sources (those 
with L x > 2 X 10 38 ergs _1 ) and the SFR of the galaxy: 

N(L > 2 x 10 38 ergs _1 ) = (2.9 ± 0.23) x SFRfMoyr" 1 ] (1) 

Using the current SFR of ~ 0.5 - l.OM yr _1 cal- 



culated by iThronson et alJ lll988l) it is predicted that in 
NGC4214 there should be between 1 and 3 sources with 
Lx > 2 X 10 38 ergs _1 . We detect one such source in this ob- 
servation, Src 11. While many ultraluminous X-ray sources 
(ULXs), which typically have X-ray luminosities of Lx ~ 
10 39 - 10 41 erg s _1 have been ob served in other st ar-forming 
galaxies (e.g. the Antennae; see lZezas et al.ll2002i) . none are 
detected in NGC 4214. 



3.3 The Luminosity Function 

Since the launch of Chandra it has been possible to study the 
X-ray point source populations in nearby galaxies in great 
detail. Studying the luminosity function of a galaxy provides 
an insight into its star formation history, and the luminos- 
ity functions of different types of galaxy can be compared 
to look for common trends. The point sources observed in 
NGC 4214 are most likely X-ray binaries. An X-ray binary 
remains active for a time which is dependent on the lifetime 
of its companion star. In the case of a high mass X-ray bi- 
nary (HMXB) the companion is an OB star with a lifetime 
of around 10 6 - 10 7 years, however the low mass compan- 
ions of low mass X-ray binaries (LMXBs) usually have a 
longer lifetime, and so LMXBs remain active f or a longe r 
period but at a lower luminosity than HMXBs JWvJl200lD . 
Luminosity functions for objects with a substantial young 
population, such as starburst galaxies, will be dominated by 
the HMXBs with no significant contribution from LMXBs, 
however as the HMXBs reach the end of their X-ray activ - 
ity the LMXBs will start to dominate (iKilgard et al.l f2002). 
Any HMXBs observed must be relatively young objects, and 
will have formed during the central starburst of NGC 4214. 
LMXBs have a lower luminosity limi t of Lx ~ 10 36 erg s _1 
<SchulJl998i . lHertz fc Grindlavlll983lL which is roughly the 
detection limit of this observation, but although some of the 
fainter X-ray sources could be LMXBs, due to the relatively 
short period since the starburst it is likely that the point 
source population is dominated by HMXBs. 

As discussed in Section [3] the predicted fluxes for the 
17 point sources within the D25 ellipse of the galaxy were 
obtained by assuming a power law model with an exponent 
of r = 1.8, assuming that they are all at the distance of 
NGC 4214. Figure shows the log(iV) - \og{L x ) plot for 
these sources. The data are best fitted with a power law 
with a slope of —0.76 ± 0.20. It can be seen that below a 
luminosity of log(Lx) ~ 36.6 the data appear to be flatter 
than the fitted slope - this is likely due to an incomplete- 
ness in the detection of the point sources due to the limits 
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Figure 5. The cumulative luminosity function for the 17 point 
sources within the D25 ellipse of NGC4214. The line depicts the 
power law fitted to the data, with a slope of —0.76 ± 0.20. 

of the instrument. The luminosity function of NGC4214 in 
comparison to those of other galaxies is discussed further in 
Section 15.11 



4 DIFFUSE EMISSION 

The faint diffuse emission from NGC 4214 can be seen in the 
centre of the Chandra image shown in Fig. Q It is clear that 
the diffuse emission comes from a region far smaller than 
the D25 ellipse, as was shown in the blow-up of the central 
region in Fig. [5] 

The X-ray contours in Fig.|2jshow that there is a bright 
ring of emission with a diameter of very roughly 0.6' (510pc). 
In addition, three point sources close to the edge of the ring 
can be seen. The diffuse emission from the galaxy can be seen 
to a greater extent in the XMM-Newton image (Fig. |3J; in 
fact the radius of this emission is around 1.4', almost twice 
that than detected by Chandra. However the spatial reso- 
lution is not good enough to clearly distinguish the point 
sources observed with Chandra from the rest of the cen- 
tral diffuse emission. Although the detection limit of the 
Chandra observation is similar to the lower luminosity limit 
of LMXBs it is possible that if several sources are present 
within a small area, for example in a cluster or associa- 
tion, they will not be distinguishable as separate individual 
sources and their emission is likely to contaminate that of 
the observed diffuse emission. 

4.1 X-ray Analysis 

The X-ray spectrum of the diffuse emission of NGC 4214 
from the Chandra observation was extracted using standard 
CIAO software, in an energy range of 0.3-8.0 keV, exclud- 
ing the three point sources from the analysis. The back- 



ground used was from the blank-sky background files asso- 
ciated with the data set - this allowed the background to be 
taken from the same position in the instrument as the source 
emission. The region used for extraction was an ellipse with 
axes of 0.64' and 0.60' (~ 550 x 500pc) chosen to encompass 
the diffuse emission. 

The diffuse emission spectrum contains 470 counts and 
was fitted with several different models in an attempt to 
obtain the best fit, each with a component, wabso, due to 
Galactic absorption fixed at Nh = 1.5 x 10 20 cm~ 2 . 

The model which provided the best fit to the data was 
wabsciwabs * mekal + wabs * mekal + po) which consists 
of two absorbed thermal components (the mekal compo- 
nents model the emission spectrum from hot gas and include 
line emission from several elements) and a power law com- 
ponent to account for any unresolved point sources. The 
statistics of the fit give a reduced % 2 °f 0.77 for 33 de- 
grees of freedom. The two thermal components both have 
absorptions of Nh = 1.9 x 10 21 cm -2 and temperatures of 
0.14 keV and 0.52 keV respectively, with the power law com- 
ponent having a photon index of T = 1.56. The unabsorbed 
flux of the diffuse emission is 3.4 x 10~ 13 erg s -1 cm -2 (0.3- 
8.0 keV), which corresponds to an unabsorbed luminosity 
of Lx = 3.5 ±3 2 o 5 xl0 38 erg s _1 (non-absorption corrected 
L x = 1.1 x 10 38 erg s" 1 ). Of this emission around 70±H% 
was due to the softer component, 16±fi% to the harder 
and 14±i % to the power law component. The errors on 
these percentage contributions show that while there is a 
degree of uncertainty as to the percentage contribution of 
the softer and harder thermal components, the division of 
the emission between the combination of the two thermal 
components and the power law component appears to be 
fairly well constrained. 

The parameters of the fit, together with the errors at 
the 90% confidence level (1.64<r) are shown in Table 01 It is 
noted that if the abundances are allowed to vary (as opposed 
to being frozen at solar) unphysical values from the fit are 
obtained - this is most likely due to there being too many 
free parameters for values to be constrained. 

The spectral response of the ACIS CCDs is known to 
vary over time due to both the evolution of the electronic 
gain of the CCDs, and the variance of their charge transfer 
inefficiency (CTI). A correction for the time dependence of 
the gain of the ACIS chips using the corr Jtgain tool 4 , was 
applied to the data. When the data were reanalysed with 
this correction it was seen that the values of the absorption 
and temperature components changed very little, although 
the slope of the power law component of the spectrum fitted 
to the diffuse emission became flatter. The effect of this cor- 
rection was small, and corrjtgain is not yet an official CXC 
product, so as such the results presented here and used for 
further analysis are those from the uncorrected data. 

Concerning the XMM-Newton data, although XMM- 
Newton is more sensitive than Chandra it is hampered by 
the fact that its spatial resolution is poorer. From the Chan- 
dra analysis it is known that there are three point sources 
(Srcs 10, 11 and 14) amongst the diffuse emission in the 
centre of the galaxy, and these could be excluded from 



contributed CXC software, see 

://asc.har vard.edu/cont-soft /software/corr_tgain.l.0.html| 
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Figure 6. The spectrum of the diffuse Chandra X-ray emission 
from NGC4214 (0.3-8.0 keV), fit with a 2 temperature mekal 
model, giving temperatures of 0. 14 keV and 0.52 keV, and a power 
law component with a photon index of V = 1.56. The contribu- 
tions of the individual components are depicted by the dashed 
lines. 




0.5 1 Z 5 



Figure 7. The spectrum of the diffuse XMM-Newton X-ray emis- 
sion from NGC 4214 (from the two MOS cameras only). The spec- 
tra was fitted over an energy range of 0.3— 8.0keV with a 2 tem- 
perature mekal model with to power law components to take into 
account unresolved point sources, one fixed with the properties 
obtained from the Chandra fits to Srcs 10, 11 and 14 (see text for 
further explanation). The softer thermal component is fit with a 
temperature of 0.23 keV, the harder with 0.59 keV. 

the Chandra analysis. However, in the XMM-Newton ob- 
servation these cannot be distinguished from the diffuse 
emission with sufficient accuracy to allow them to be ex- 
cluded without removing significant parts of the diffuse emis- 
sion as well. To take into account these 3 sources an addi- 
tional absorbed power law component with values frozen at 
Nh = 1-8 x 10 21 cm~ 2 and T — 1.8 was included in the mod- 
els (the values used were based on those obtained from the 
fits to Srcs 10 and 11 - see Table 0. The spectra were again 
fitted with a wabso(wabs * mekal + wabs * mekal) model, 
with the additional component to take into account the cen- 
tral point sources included. It was found that the spectra 
from the two MOS cameras showed good agreement with 
the results obtained from Chandra; however, when the re- 
sults from the PN were included a fit with unphysical values 
was obtained. The discrepancy between the different cam- 
eras is most probably due to the differences in calibration 



Table 4. The parameters from the best fit model to both the 
Chandra and XMM-Newton data. Each data set was fitted with a 
wabsQ * (wabs *mekal+wabs *mekal+po) model; the XMM-Newton 
data had an additional power law component with an exponent 
of T = 1.8 and Njj = 1.8 X 10 21 cm -2 to take into account the 
unresolved central point sources (Srcs 10, 11 and 14) - see text. In 
each case wabsQ was assumed to be the Galactic absorption value 
of Nh = 1-5 X 10 20 cm~ 2 . The errors shown are those calculated 
at the 90% confidence level. 



Chandra Data XMM-Newton Data 



N H i (10 22 cm- 2 ) 


n i q_i_0.53 
u ' ia=t 0.19 


n 17 _i_0.06 
U - L ' :t 0.03 


A;Ti (keV) 


n 14+0. 06 


0-23±0.04 


N H 2 (10 22 cm- 2 ) 


n in+0.48 
u ' iy=t 0.19 


0.85±g : 2 3 7 


kT 2 (keV) 


n w+0 -19 
U.OZrtg 44 


o.59±g:i§ 


r 


1.56±?If 


1 OQ + 0.51 


X 2 /d.o.f. 


0.77/31 


1.11/40 


L x (lO^ergs- 1 ) 


3.53±3 2 5 


k ru_i_l.37 
D - U * =c 1.64 



between them. The fits were thus carried out using the MOS 
cameras only. 

The fit to the XMM-Newton data gives two absorp- 
tion components of Nh = 1.7 X 10 21 cm -2 and Nh = 
8.5 x 10 21 cm~ 2 , a temperature for the softer component of 
kT = 0.23 keV and for the harder of kT = 0.59 keV and 
a power law component with a photon index of T = 1.93. 
The results of the fit, together with the errors on the fitted 
values are given in Table |1] The unabsorbed luminosity was 
found to be Lx = 6.0 ±};| xl0 38 erg s -1 , higher than that 
obtained from the Chandra fit which is to be expected due 
to the inclusion of point source emission. Indeed, when the 
calculated contribution of Srcs 10, 11 and 14 is subtracted 
from this value a luminosity of Lx ~ 2.5 x 10 38 erg s _1 is 
obtained which is in better agreement with that from the 
Chandra observation. The contribution of the various com- 
ponents to this emission was found to be 19 ± 2 % for the cold 
component, 34 ±11 % for the hot, 9 if % due to the power 
law component, and 38 ±3 % due to the Chandra sources 10, 
11 an 14. The results from the XMM-Newton data appear 
to broadly confirm those obtained from Chandra. 

By making some basic assumptions it is possible to 
calculate additional parameters of the gas components of 
NGC 4214. The following calculations were carried out 
using the results obtained from the Chandra analysis. 
iMafz-Apellaniz et alJ (Il999l) suggest that NGC 4214, which 
we are observing face on, has an unusually thin disk with a 
total height of around 200pc extended by outflows by per- 
haps a further 200pc. As such a cylindrical geometry for the 
galaxy is assumed with a thickness of 200pc and a diameter 
of ~ 0.6' (510pc) which corresponds to that of the observed 
emitting region. 

The mass of the emitting gas can be calculated from the 
normalisations of the model components. The normalisation, 
k is given by 

10~ 14 f 

where D = 2.94Mpc, V is the volume of the emitting region, 
and n e and n H are the number densities of electrons and 
hydrogen ions respectively. We assume that the gas is fully 
ionised (n c ~ n n)- The total mass of the emitting gas is 
~ 2.1 x 1O 5 77 1//2 M0 with the two thermal components and 
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Table 5. Parameters calculated from the best fit properties from 
the Chandra analysis of the diffuse emission. For assumptions 
used in the calculations see text. 



Emission 


Softer 


Harder 




component 


component 


kT (keV) 


0.14 


0.52 


T (K) 


1.62x10 s 


6.04x10 s 


L x (ergs 1 ) 


2.48X10 38 


5.57X10 37 


n e (cm- 3 ) X (1/y/fj) 


0.10 


0.03 


E Th (erg) x{y/rj 


1.44x10 s3 


1.61x10 s3 


M (M S ) x (v^T) 


1.82x10 s 


5.50X10 4 


P (dyn cm -2 ) X (1/y/rj) 


4.61 x 10" 11 


5.16 x 10" 11 


tcooi (yr)x(y^) 


1.86 XlO 7 


9.32X10 7 



one power law component contributing 1.5 x IO^^Mq, 
4.2 x 1OV /2 M and 1.8 x lO 4 r? 1/2 M respectively. We note 
the dependence on rj, th e volume filling factor of the X-ray 
emitting material - see IStrickland fc Stevens! i2000l) for a 
discussion of its pro perties. We make no assum ptions about 
th value of n, but IStrickland fc Stevens! ll200fj| find that it 
could be quite low (~ 0.1) for soft emission; if this is the case 
both the mass and energy content of the gas are reduced. 

Further parameters of the emitting gas can be calcu- 
lated by making some basic assumptions about its prop- 
erties. Using the emission integral, EI, from Equation 
((fc x 47rD 2 )/10~ 14 ), the electron density can be calculated 
using n e ~ (EI /Vt)) 1 ^ 2 , assuming the emitting volume, 
V, calculated previously. The dynamical pressure, P, of 
the gas is obtained using P ~ 2n e kTn~ 1 ^ 2 , and the ther- 
mal energy, Eth, of the gas from Em ~ 3n e kTVrj ll/2 . Fi- 
nally the cooling time, t coo i, of the gas is calculated using 
tcooi ~ 3kTn 1/2 /An e where A = L X /EI. The results of 
these calculations for the two components of the gas are 
shown in Table |S] 

4.2 Comparison with Optical Observations 

As has been noted before, excluding the point sources, 
the bulk of the X-ray emission from NGC4214 comes 
from the central regions of the galaxy. The centroid of 
the X-ray emission (located at a[2000] = 12 h 15 m 39.77 s , 
5[2000] = 36°19'28.9", see Fig. coincides with the NW 
star-form ation complex, marked out by Ha emission, ob- 
served bv lMacKentv et alJ fcOOfl) . 

Fig . El shows th e Ha emission from NGC4214 observed 
bv lMacKentv et al.1 (2000) with the X-ray contours from the 
Chandra data overlaid. This image shows clearly that there 
is more X-ray emission associated with the NW star-forming 
complex than with the SE complex. This is perhaps to be 
expected as the NW complex is the more evolved and more 
massive of the two, hence it contains more supernovae, and 
any superbubble will have had more time to develop. Al- 
though the bulk of the X-ray emission is centred around 
the point sources in the NW complex there is emission ex- 
tending away in a southerly direction which gives the im- 
pression of having a ring-like shape. This emission is most 
likely solely due to the ongoing star formation in the NW 
complex, with the shocked gas from stellar winds and su- 
pernovae preferentially expanding into regions of the ISM 
with lower ambient density. The geometry of the emission 



is somewha t reminiscen t to that of a so-called 'champagne 
phase' (see lYorkelll986l for a review). There is little X-ray 
emission corresponding with the position of the SE complex, 
which is further conformation that this is the younger of the 
two star formation regions. It appears that substantial su- 
perbubbles have not yet had time to form, and hence the 
region is likely not yet able to inject significant energy into 
the surrounding medium. 

Src 10 corresponds to the position of a star cluster be- 
lieved to be the dynamical centre of the galaxy (as dis- 
cussed in more detail in Section 13.21 . A young SSC, of 
age ~ 3.0 — 3.5 Myr, has been identified in optical im- 
ages at a [20001 = 12 /t 15 m 39.44 s , (5120001 = 36°19'35.0" 
llieitherer et alJl996tlMacKentv et al . 2000), and it appears 
to reside in an Ha cavity. At this position there is increased 
emission in the X-ray suggesting the presence of cluster type 
emission. 



5 DISCUSSION 
5.1 Point Sources 

The luminosity functions of galaxies can now be studied in 
great detail due to the fact that Chandra is able to detect 
many individual X-ray point sources within the galaxies. 
The slope of the luminosity function can be used as an in- 
dicator of some of the general properties of the host galaxy. 
A recent study o f 32 n earby spiral and elliptical galaxies 
bv lColbert et alJ (|2003) shows a clear distinction between 
the slopes of the luminosity functions of these two types of 
galaxies, with those of the spiral galaxies being much steeper 
than those of the ellipticals. 

In those galaxies where a large number of point 
sources has been observed the luminosity functions arc 
often fitted with po wer laws, f or example NGC 444 9 
JSummers et al]|2003al) NGC 5253 dSummers et al.||2003rj|) . 
M31 JSupper et alJll997h . and M81 jTennant et al.ll200ll) . 

If a continuous uniform star formation process occurs 
then the luminosity function of the sources should appear as 
unbroken, but when a starburst is triggered new HMXBs will 
be formed which will appear as a break in the function. It is 
therefore not unlikely that a break in the luminosity function 
will be present in all starburst galaxies, the position of which 
will decrease in luminosity with time, perhaps providing an 
indicator to the time of previous bursts in the galaxy. It is 
possible that no break is seen in NGC 4214 due to the limited 
statistics and the small number of point sources observed; if 
there is contamination of the luminosity function by LMXBs 
these could mask a break in the HMXB population when 
there is such a small sample being considered. 

The results for the X-ray luminosity function presented 
here for NGC 4214 can be combined with those of other 
dwarf starbursts and then further compared with results 
from a wider sample of galaxies for a range of parameters 
(for example iKilgard et alJ |2002L see below). The relevant 
parameters (and sources) for each galaxy used in this study 
are shown in Table 

The blue luminosity Lb was calculated using iTullvl 
(1988) 

logL s (L Q ) = 12.192 -0AB T + 2\ogD (3) 
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Figure 8. The Chandra observation of the central regions of NGC4214. X-ray contours from the binned, smoothed image are super- 
imposed onto an Ha image. In places regions of enhanced X-ray emission clearly correspond to regions of increased Hq emission, but 
that in others there appears to be no correlation. This H« image was kindly provided by Jesus Mai'z-Apellaniz from data published in 
MacKenty et al. (2000) and also available from http://www.stsci.edu/~jmaiz 



with D in Mp c, and Lfir calculated using the 60fj,m and 
100/xm fluxes dPevereux fc Ealedll989F) 

L FIR {L e ) = 3.65 x 10 5 [2.58S 6 o M ™ + S l0 ^ m ]D 2 (4) 

iKilgard et al] J2002) examined the point source lumi- 
nosity functions of 7 galaxies; 4 spirals and 3 starbursts. 
The 3 'normal' starbursts, M82, NGC 253 and the Anten- 
nae (NGC 4038/9), have flatter luminosity functions than 
the spirals, with slopes of —0.50, —0.81 and —0.53 respec- 
tively compared to a range of —1.07 to —1.30 for the spi- 
rals (see Table |5J. We supplement this sam ple with results 
for the dwarf starbu rsts galaxies NGC 4449 dSummers et al] 
l2003al) . NGC 5253 iSummers et al.ll2003bJl and NGC 4214 
(this paper). 

The far infrared luminosity Lfir of a galaxy can be 
used as an indication of a young stellar populati on and of 
its current rate of star formation ]Kennicuttlll99gf) . Thus, if 
the luminosity function slope is dependent on the rate of star 
formation, then we might expect a correlation between the 
luminosity function slope and Lfir- Fig. |U] shows the lumi- 
nosity function slope plotted versus Lfir, and it is obvious 
that the starburst galaxies have a flatter slope than the 4 
normal spiral galaxies, suggesting that the higher the star 
formation rate the flatter the luminosity distribution. This is 
probably due to the fact that there is a larger proportion of 



young bright sources (e.g. HMXBs ) in star-b ursting regions 
than is present in spiral galaxies (IWull200lD . However, the 
dwarf starbursts also have lower luminosity function slopes, 
but because the dwarf starburst galaxies have Lfir values 
comparable to the spirals there is no clear segregation on 
the basis of Lfir- 

A clearer segregation can be seen in Fig. 1101 which 
shows the X-ray luminosity function slope plotted against 
the ratio of the 60/xm (Seo/im) to lOO^m (Siooum) flu x for the 
spiral and starburst galaxies. iLehnert fc Heckmanl (^996a) 
state that if Saoum/ Siooum > 0.5 then it is more likely that 
a galactic superwind will occur. As is clear from Fig. 1101 
this criterion is satisfied for all the starburst galaxies in the 
sample. A superwind is indeed observed in both M82 (e.g. 
IStrickland etail Il997l ILehnert et all Il999h and NGC 253 
(e.g lStrickland et al.ll2000D . In NGC 5253 multiple super- 
bubb les are observed in the X-ray iStrickland fc Stevensl 
1999) and Ha filament s are seen extending away from the 
galaxy jMarlowe et al1ll995l) . F rom NGC 4449 both X-ray 
JSummers et al]l2003al) and Ha jMarlowe et aljll995l) emis- 
sion from a superbubble is observed, but a large HI halo 
surrounding the galaxy may prevent the escape of mate- 
rial into the surrounding intergalactic medium. Using the 
Chandra X-ray data it has been calculated that NGC 4214 
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is likely to undergo a blow-out; this is discussed in more de- 
tail in Section lo.3l As the Antennae are two merging galaxies 
as yet no superwind has formed due to the disruption from 
the ongoing interaction. 

Fig, II II shows the luminosity function slope against the 
60^im flux scaled to the area of the D25 ellipse as a mea- 
sure of the star formation rate per unit area. Unsurpris- 
ingly, this tends to be higher for the starburst galaxies than 
for the spirals (with the exception of M83), and again we 
see a segregation in luminosity function slope with star- 
forming activity. M83 has been classified as a nuclear star- 
burst galaxy, and HST observations show a lar ge number of 
youn g, massive clusters in the nuclear region dHarris et alJ 
l200ll) . lKilgard et al.l ( I2002T) state that most of the detected 
point sources used in the analysis of the luminosity function 
of M83 are from regions outside the central star-bursting 
region. The fact that only the disk population is analysed 
means that it can be grouped with the other spirals in 
the sample when considering the luminosity function slope. 
However, other values such as the optical luminosity (Lb) 
and the far infrared luminosity (Lfir) will have this central 
region included which will affect the results. Indeed, inspec- 
tion of the 60/im IRAS image of M83 shows a large variation 
of emission across the galaxy, with the highest contribution 
coming from the centre. For analysis purposes, however, the 
flux from the entire galaxy was used. This explains why M83 
seems to agree with the properties of spirals in respect to 
the luminosity function, but for other properties relating to 
the whole galaxy it is more similar to what is seen in the 
starbursts. The luminosity function of M82, another star- 
burst galaxy with localised star formation, contains point 
sources from the entire galaxy including the regions of star 
formation. Unlike M83 the luminosity function slope of M82 
shows agreement with those of the other starburst galaxies. 

Typically Lfir scales with Lb for normal galaxies, but 
in strong starbursts Lfir is seen to increase more rapidly 
than Lb- Fig. 1121 shows the slope of the luminosity func- 
tion plotted against \og[LFiR/ Lb]- In general, the starburst 
galaxies in the sample have a higher [Lfir/ Lb] ratio, which 
is to be expected, but if we exclude M83 again there is still 
evidence of a change in luminosity function slopes at a value 
of log Lfir/ Lb ~ —0.75. M83 has the highest [Lfir/ Lb] 
ratio of the spiral sample; this could again be due to the fact 
that it has had some recent star formation. 

We note in passing that we can compare the total num- 
ber of sources in each galaxy with log Lx > 37.5 (a lumi- 
nosity level which provides a compromise in that sources 
this faint are detected in more distant and massive galax- 
ies such as the Antennae and sources this bright are de- 
tected in smaller and nearby dwarfs such as NGC 4214 and 
NGC5253). If we plot this number against Lb then we find 
a linear correlation with on average 0.2 — 1 X-ray sources 
with L x > 10 37 ' 5 erg s _1 for each 10 8 L© in Lb- 

It is worth noting that NGC 4214 was analysed with a 
correction for the low energy quantum efficiency degrada- 
tion of Chandra (see Section l2.ll . which gave slightly dif- 
ferent values than when it was analysed without this cor- 
rection. The rest of the data was obtained from previously 
published material which was analysed before the correction 
was available and it is not known how much this will have 
affected the validity of the results. However, this problem 
is getting worse with time and will therefore have had less 
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Figure 9. The luminosity function slope plotted against Lfir 
for the lKilgard et all 120021) sample of 4 spirals and 3 starbursts, 
with the results for 3 dwarf starbursts (NGC 4214, NGC 4449 and 
NGC 5253) added. The spiral galaxies are denoted as circles in the 
plot; the starbursts as stars. 
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Figure 10. The luminosity function slope plotted against the 
ratio of the 60/xm to 100/xm flux for 4 spirals (circles) and 6 star- 
bursts (stars). A division can be seen between the two types of 
galaxies, both in the luminosity function slope and in the flux 
ratio. 



of an effect on observations such as these which were taken 
early in Chandra's lifetime, and so it is assumed that it is 
still valid to compare the results obtained. 
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Table 6. Luminosity function data for a sample of nearby galaxies; 3 dwarf starbursts, 4 starbursts and 3 spirals. Much of the data 
is from the sample of lKilgard et alj (2002), along with the data from NGC4214, NGC4449 and NGC5253 mentioned in the text. L B 
and Lpjn are estimated using Eqn. 131 and Kl The galaxy area is estimated from (n ab)cos8 where a, b are the galaxy semi-major and 
semi-minor axes of the D25 ellipse, and 8 is the inclination of the galaxy, obtained from the LEDA database. 





Distance 


LI slope 


oU/im iiux 


lUU/im nux 




logL FIR 


logL B 


-D25 area 


Galaxy 


(Mpc) 




(Jy) 


(Jy) 




(erg s _1 ) 


(erg s _1 ) 


(arcmin 2 ) 


NGC 4214 


2.9 


0.70 


17. 9 a 


29. a 


0.62 


41.96 


42.62 


29.2 


NGC 4449 


2.9 


0.51 


36.0 6 


73.0 6 


0.49 


42.30 


42.72 


5.3 


NGC 5253 


3.2 


0.54 


10.5 C 


29.4 C 


1.04 


41.90 


42.43 


3.3 


M82 


5.2 


0.50 


1313.5 a 


1355.4 a 


0.97 


44.26 


43.49 


4.0 


NGC 253 


3.0 


0.81 


931. 7 a 


1861. 7 d 


0.50 


43.73 


43.52 


31.1 


Antennae 


25.5 


0.51 


48.7° 


82.0" 


0.59 


44.28 


44.15 


12.6 


NGC 1291 


8.6 


1.07 


1.8 d 


10.1 d 


0.17 


42.18 


43.89 


70.3 


M83 


4.7 


1.38 


266.0 d 


638.6 d 


0.42 


43.61 


43.84 


110.8 


NGC 3184 


8.7 


1.11 


8.9 a 


29.0 a 


0.31 


42.74 


43.51 


32.0 


IC 5332 


8.4 


1.30 


0.8 C 


5.1 c 


0.16 


41.85 


43.12 


25.3 
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Figure 11. The luminosity function slope plotted against the 
60/jm flux divided by the area of the D25 ellipse. (Spiral galaxies 
are denoted as circles; starbursts as stars). 



5.2 Star Formation Rate of NGC 4214 

iRanalli. Comastri fc SettH ll200ot) propose that the X-ray lu- 
minosities of star-forming galaxies can be used as an indica- 
tion of the SFR of the galaxy. They derive the following re- 
lation between the 2-10 keV X-ray luminosity and the SFR: 

SFR(M B yr- 1 ) = 2.0 x lO^La-iokcv (5) 

The spectrum of the total X-ray emission within the 
-D25 ellipse of the galaxy was extracted from the Chan- 
dra data and fitted with an absorbed thermal component 
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Figure 12. The luminosity function slope for 4 spirals (cir- 
cles) and 6 starbursts (stars) plotted against \og(Lpm/ Lb)- This 
shows that in this sample although in general the starbursts have 
a higher Lpju/Lb ratio there is no clear division between the 
two types of galaxies. 



plus power law model within XSPEC. The absorbing col- 
umn density was divided into 2 components, with one frozen 
at the Galactic column density Nh = 1.5 x 10 20 cm -2 and 
the other allowed to fit. The exponent of the power law 
component was frozen at a value of T = 1.8 to take into 
account the point sources emission. The fitted column for 
the additional absorbing component, local to NGC 4214, 
was found to be Nh = 1.4 x 10 21 cm -2 , giving a total col- 
umn of 1.6 x 10 21 cm -2 . This basic spectral analysis pro- 
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vided an unabsorbed 2-10 keV X-ray luminosity ol Lx 



3.4 x 10 38 ergs 



When combined with the relation from 

for 



iRanalli et all (2003) this gives a value of 0.0 7Mq yr 
the SFR of NGC 421 4. 

iKennicuttl il998fl derives the following relation between 
SFR and the FIR luminosity in a galaxy: 



SFR(M Q yr- 1 ) = 



2.2 x 10 43 



(6) 



Using a value of Lpm = 9.1 x 10 41 ergs _1 (estimated using 
Eqn.^J the Kennicutt relation gives a SFR of O.O4M yr" 1 , 
in good agreement with the value obtained using the X-ray 
properties of the galaxy. 



5.3 Diffuse X-ray Emission and IR Properties 

In this observation we have been able to resolve point sources 
and diffuse emission in the central regions of NGC 4214. This 
in turn enables us to reliably estimate, for the first time, the 
energy present in the ISM and ultimately whether we expect 
an outflow from NGC 4214. 

Metals and gas are potentially ejected from dwarf star- 
burst galaxies into the surrounding intergalactic medium, 
powered by the strong stellar winds from OB stars, and 
then after ~ 3.5Myr predominantly by the energy injected 
by supernovae, which adds to the mass loss from the g alaxy. 
Indeed, observations of NGC 1569 ijMartin et al]E o02'l show 
that the diffuse X-ray emitting gas has a supersolar a/Fe ra- 
tio. These results suggest that the hot diffuse gas has been 
enriched by the ejecta of explosions Type II supernovae in- 
jecting heavy elements into the surrounding medium. The 
winds from the stars and supernovae in clusters and stel- 
lar associations will combine, sweeping up the surrounding 
interstellar medium, creating a superbubble of hot, X-ray 
emitting gas. Surrounding the superbubble will be a cooler 
shell of material which emits Ha line emission. If instabilities 
within the bubble occur they can cause the shell to rupture, 
allo wing gas to escape in the form of a superwind. 

IWeaver et alJ 1)19771) model bubbles for single stars, and 
these results are often used to model the formation of bub- 
bles around clusters and associations that contain many 
massive stars. However, this model must be used with cau- 
tion. In the case of NGC 4214 there are several separate 
superbubbles that have been observed, which correspond to 
the dif ferent clusters and asso ciations at the centre of the 
galaxy. iMacKentv et all J200CT) observe bubbles around the 
star forming regions in the NW com plex that are a round 
30% smaller than expected using the IWeaver et aD lll977l) 
model, an d conclude that the model cannot be applied in 
this case. Summers. Stevens fc Strickland! (|200lfr find that 
the predicted expansion velocity for the dwarf starburst 
galaxy Markarian 33 from using the Weaver model is less 
than that observed. 

In their observations of NGC 4214 lMafz-Apellaniz et aD 
(1999) observe no bubbles surrounding the three star clus- 
ters in the SE complex. They suggest that observations 
of the SE complex do not agree with a blow-out having 
occurred, and that despite the Weaver model predicting 
that there should be a superbubble with Rb ~70pc and 
vb ~14km s~ it is likely that a su perbubble does not ac- 
tually exist. IMacKentv et alJ i2000t) suggest that there is a 
delay in the formation of a superbubble in clusters of about 



2 Myr due to the time it takes for the high number of individ- 
ual stellar bubbles produced to combine and break out of the 
surrounding medium. If this is indeed the case it would ex- 
plain why X-ray emission is seen to predominantly coincide 
with the position of the somewhat older NW complex, where 
cluster superbubbles have had time to form, but not with 
the SE complex where there is a lack of shocked, hot gas. 
The time delay in the formation of these clusters could have 
an impact on e stimations of the amount of energy released 
by the galaxy (M afz-Apellaniz et alj fl999). The mechani- 
cal energy could be dissipated via the heating of interstellar 
gas, which would lead to an overestimation of the energy 
released, as would the fact that stellar synthesis models as- 
sume that the massive stars would appear earlier than is in 
fact the case. 

iMac Low fe McCravl (ll988T> derive the following criteria 
for whether a blow-out from a galaxy is to be expected: 



A-, 



Jmech,i\ H kpc P 4 



-3/2 1/2 



(7) 



where L mec h,4i is the mechanical energy luminosity in units 
of 10 41 erg s _1 , Hkpc is the scale height of the galaxy 
in kpc, P4 is the initial pressure of the ISM in units of 
P/k = 10 4 K cm and no is the initial density of the 
ISM. A mm is a measure of the rate of kinetic energy in- 
jection, and if A mm > 100 then a blow-out will occur. From 
Table the total thermal content of the X-ray emitting 
gas is around 3.1xl0 53 erg s _1 . Assuming that the observed 
diffuse X-ray emission comes from the most recent burst 
of star- formation ass ociated with the NW complex, with 
an age of ~ 3.5Myr l|MacKentv et al1l200ffl . we estimate 
that L mec h ~ 2.8 x 10 39 erg s _1 . The initial pressure of the 
ISM is assumed to be equal to that of the Milky Way, i.e. 
Pi ~ 1. For plausible values of the initial gas density no 
(such as the initial density being 10 x greater than the cur - 
rent derived densities in Table||l cf ISummers et alJl2003^ 
for Amm > 100 and a blow-out to occur then we require 
Hkpc ;$ 3kpc. Obviously the derived upper limit for the scale 
height necessary for a blow-out to occur is very weakly de- 
pendent on the assume d value of no. The assumed pressure 
is comparable to that of lMac Low fc Ferraral lll999l) who as- 
sume Pj ~ 4 x 10 3 Kcm~ 3 . Indeed, if the lower value of 
IMac Low fc Ferraral l)l999Fl is used it makes little difference 
to the result (Hk pc < 2kpc). The estimated value of Hk pc is 
significantly large r than the optically observed dimensions of 
the galactic disk (|Mai'z-Apellaniz et al.lll999l) . and as such 
we conclude that in this galaxy a blow-out from the disk 
ISM is likely to o ccur. 

We note that lWilcots fc Thurowl (|200lfl calculate Lmech 
to be ~ 10 39 ergs _1 , a factor of ten smaller than the value 
obtained from our analysis. However, if Lmech is indeed a 
factor of ten smaller this will reduce the scale height by a 
factor of \/T0, which will still be larger than the observed 
size of the optical disk. 

iLehnert fc Heckmanl jl996al) have studied a large sam- 
ple of starburst galaxies, looking at which properties make 
them conducive to galactic superwinds being formed, and 
find that those galaxies with extreme IR properties are far 
more likely to have a superwind present. One of the crite- 
ria is that the galaxy has a ratio of Sea^m/ Sioa^m. > 0.5; 
NGC 4214 has a value of 0.62, implying that a superwind 
is a possibility. However, NGC 4214 fails to meet the other 
two criteria; that Lfir > 10 44 ergs _1 and Lfir/Lb > 2. 
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These results appear to suggest that there is not a high 
chance of this galaxy forming and supporting a galactic su- 
pe rwind, however it should be noted that the galaxies used 
bv lLehnert fc Heckmanl (ll996aT) to derive these criteria were 
mainly spiral galaxies, and as such the results may not be 
quite so relevant for a dwarf starburst such as NGC 4214. 

In the sample of 9 galaxies used in this to study the 
properties of luminosity functions of starbursts, all of the 
starbursts have Ssoum / 'Sioofim > 0.5 (note that NGC 4449 
has a value of 0.49). However, with the exception of M82, 
with Lfir/Lb = 5.9, none meet the other criteria for the 
formation of a superwind to be likely. 



6 SUMMARY AND CONCLUSIONS 

We have presented the results from both Chandra and 
XMM-Newton observations of NGC 4214, which show the 
diffuse X-ray emission from this dwarf starburst galaxy. The 
high spatial resolution of Chandra has allowed the detailed 
structure of the X-ray emission of NGC 4214 to be observed 
for the first time. 17 point sources were observed within 
the D25 ellipse of the galaxy, and further analysis of these 
suggest that the majority are HMXBs associated with the 
galaxy as opposed to being background AGN. 

The luminosity function of the point sources is simi- 
lar to those seen in other starbursts such as NGC 5253 and 
NGC 4449. Indeed, when combined with thos e of NGC 4449 
and N GC 5253 and added to the sample of iKilgard et alJ 
(2002) this data strengthens the argument that the lumi- 
nosity functions of starburst galaxies are flatter than those 
seen in normal spirals, and indeed extends this result to in- 
clude dwarf starburst galaxies. It is seen that the starburst 
galaxies tend to have a higher 60[im to lOO/xm flux ratio than 
spirals, a cr iteria for being more likely t o produce a galactic 
superwind jLehnert fc Heckmanl fl996aTl . However, there is 
no clear segregation between the spirals and the starbursts 
when the luminosit y function slopes ar e comp ared to the 
ratio of Lftb/I/h. iLehnert fc Heckmanl lll996al) state that 
for a superwind to form a large excess of IR luminosities is 
needed, and although this criteria is not met in our sample 
overall the starbursts show a larger Lfir/Lb ratio than the 
spirals, as would be expected. In general the starbursts have 
a larger amount of star formation per unit area (using the 
D25 ellipse area as a scaling parameter). 

The diffuse emission from the galaxy is observed in both 
the Chandra and XMM-Newton observation, and is seen to 
almost twice the extension in XMM-Newton (radius ~ 1.4') 
than in Chandra. The data suggest the presence of gas at 
two separate temperatures, with Chandra analysis provid- 
ing values for these of around 0.3 and 4.6 keV and a com- 
bined absorption corrected flux of 3.4 x 10~ 13 erg s~ 1 cm -2 
corresponding to an unabsorbed luminosity of Lx = 3.5 x 
10 38 erg s _1 . Using derived X-ray properties the SFR of 
NGC 4214 is calculated as O.O7M0 yr _1 , in reasonable agree- 
ment with the value of O.O4M0 yr _1 derived using Lfir- 

The X-ray emission coincides with regions of Ha emis- 
sion and hence is coincident with star forming regions. The 
presence of X-ray emission coincident with the NW complex, 
but not the SE complex of Ha emission agrees with previous 
speculation that there is a time delay before the wind from a 
cluster of stars can break out and start interacting with the 



surrounding medium. Using some basic assumptions about 
the properties of the X-ray emitting gas is has been calcu- 
lated that it is likely that a blow-out will occur from the 
disk of NGC 4214, enriching the surrounding medium with 
materials from this starburst. 
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